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Abstract 


Precipitation is one of the main causes of floods which due to the improper dispersion of 
rain gauge stations and the newly established ones in most basins of the country, the use 
of the precipitation data faces serious challenges. Remote-sensing method can be used in 
such cases. In this study, the reanalysis rainfall data of ERAS for Kashafrood catchment 
were evaluated in daily and monthly timescales and then its streamflow hydrograph was 
evaluated using rainfall data of Zoshk station and Zoshk watershed parameters in the 
HEC-HMS software. The reanalysis rainfall data of ERAS showed underestimation and 
overestimation and also had better accuracy in monthly timescale (R = [0.87-0.99]) than 
daily timescale (R = [0.89-0.62]) and its flow hydrograph compared to the observational 
hydrograph was more satisfactory (NSE = 0.49 and PBias = -26.46). Flood zone of 
reanalysis precipitation data of ERAS in relation to the reference flood zone was -0.07, 
which was due to the underestimation of precipitation data and the subsequent 
underestimation of its hydrograph peak flow in the HEC-HMS model. Future studies can 
be used to evaluate the flood zone resulting from precipitation data of various products in 
HEC-RAS model. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. شمارة چهل و سوم پاییز ۶۰۱ صص ۲۳-4۲ 


مقاله پژوهشی 


صحت سنجى يهنه و هیدروگراف سیلاب با استفاده از داده‌های بارش ERAS‏ (مطالعه موردی: o>‏ ضه زشک) 


سعیده پورانتظاری - دانشجوی کارشناسی ارشدء گروه علوم و مهندسی آب» دانشگاه فردوسى مشهد. مشهد. ايران 
sD‏ اسماعیلی '- دانشیار گروه علوم و مهندسی آبء دانشگاه فردوسی مشهد. مشهد. ایران 
علیرضا فریدحسینی- دانشیار گروه علوم و مهندسی آب. دانشگاه فردوسی مشهد. مشهد. ایران 
الهه غفاری - دانش‌آموخته دکتری گروه علوم و مهندسی آب. دانشگاه فردوسی مشهد. مشهد. ايران 


تاريخ دریافت: ۱٤٠١/۱۰/٤‏ تاريخ بازنگری: ۱2۰۰/۱۱/۲۶ تاريخ تصویب: ۱2۰۰/۱۱/۲۹ 


حکیده 


بارش یکی از اصلی‌ترین عوامل وقوع سیلاب‌هاست که به علت عدم پراکندگی مناسب ایستگاه‌های 
باران‌سنجی و تازهتأسيس بودن برخی از اين ایستگاه‌ها در اکثر حوضه‌های کشور استفاده از داده‌های بارش 
با چالشی جدی روبروست. ازاین‌رو روش‌های دورسنجی می‌توانند یکی از گزینه‌های مورداستفاده در اين 
زمینه باشند. در پژوهش حاضر داده‌های بارش بازتحلیل‌شده ERAS‏ برای حوضه آبریز کشف‌رود در 
گام‌های زمانی روزانه و ماهانه مورد ارزیابی قرار كرفت و سپس هیدروگراف و پهنه سیلاب آنء با استفاده 
از داده‌های بارش ایستگاه زشک و پارامترهای مربوط به حوضه آبخیز زشک در نرم‌افزارهای HEC-‏ 
5 و HEC-RAS‏ موردبررسی قرار گرفت. نتایج نشان داد که داده‌های بارش 118۸5 دارای 
کم‌برآوردی می‌باشد و داده‌های بارش در گام زمانی ماهانه ([۰,1۸-۰,۹۲]-) دقت بالاتری نسبت به 
روزانه ([۰,0۲-۰,۰۸۱]-16) داشت و هیدروگراف جریان آن در مقایسه با هیدروگراف مشاهداتی 


غیرقابل قبول NSE=- tY)‏ و (PBias=-00,11‏ بود. ag‏ سیلاب دبی داده‌های بارش بازتحلیل شده 
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Vt‏ جغرافیا و مخاطرات محیطی شمارة سوم 


ERAS‏ نیز نسبت به پهنه سیلاب دبی مرجع معادل ۰,۱۹- به دست آمد که به دلیل کم‌برآوردی داده‌های 
بارش و به دنبال Ol‏ کم‌برآورد شدن دبی اوج هیدروگراف of‏ در مدل HEC-HMS‏ اين نتيجه به دست 
آمده است. مطالعات آینده می‌توانند در حصوص ارزیابی پهنه سیلاب حاصل از داده‌های بارش محصولات 
مختلف در نرم‌افزار HEC-RAS‏ باشند. 

کلیدواژه‌ها: پهنه‌بندی سیلاب, رودخانه زک HEC-RAS HEC-HMS ERAS‏ 


۱- مقدمه 

سیلاب آسیب‌ها و خسارات فراوانی به جان و اموال مردم وارد می‌کند. در وقایع سيل ۱۳۹۸ بنابر گزارش سازمان 
مديريت بحران کشور بیش از 1۰ هزار واحد شهری و روستایی تخریب شدند و بیش از ۷۵ هزار واحد مسکونی 
شهری و روستایی آسیب دیدند. در اين حوادث. ۳۲۸۵ نفر مصدوم و بیمار به مراكز و واحدهای درمانی مراجعه 
کردند و ۷۱ نفر از هم‌وطنانمان جان خود را از دست دادند (هیئت ویژه گزارش ملی سیلاب‌های ایران, ۱۳۹۸). بارش 
یکی از اصلی‌ترین عوامل وقوع سیلاب‌های ناگهانی در حوضه‌های آبخیز کوچک است (لو" و همکاران» ۲۰۲۲) که 
به علت عدم پراکندگی مناسب ايستكاههاى باران‌سنجی و تازه تاسیس بودن برخی از اين ایستگاه‌ها در اكثر 
حوضه‌های کشور استفاده از اين داده‌ها با جالشى جدی روبروست. ازاينرو استفاده از روش‌های دورسنجی 
می‌تواند یکی از گزینه‌های مورداستفاده در مدل‌های هیدرولوژیکی (برای دستیابی به مشخصات جریان از قبیل میزان 
حجم دبی اوج و زمان رسیدن به دبی اوج جریان) و هیدرولیکی (برای دستيابى به مشخصات سیلاب از قبیل عمق و 
تراز سطح آب و سرعت جریان و درنهایت پهنه درگیر سیلاب) باشد. داده‌های بارش بازتحلیل‌شده BRAS‏ پنجمین 
نسل از داده‌های پایگاه مرکز پیش‌بینی وضعیت آب‌وهوای اروپا (ECMWE")‏ هستند که از ۱ ژانویه 1۹۷۹ کل کره 
زمين را پوشش می‌دهند و دارای توان تفکیک مکانی و زمانی به ترتیب" ۰۰,۲۵ "۰,۲۵ و ۱ ساعت می‌باشند (سباه " و 
همکاران» ۲۰۲۱). داده‌های بازتحلیل‌شده از ترکیب نتایج حاصل از پیش‌بینی‌های کوتاه‌مدت مدل‌های عددی 
پیش‌بینی وضع هوا با انواع داده‌های محاسباتی و مشاهداتی از سنجنده‌های مختلف به‌دست می‌آیند (کریمی و 
همکاران, ۱۶۰۰). HEC-HMS‏ یکی از نرم‌افزارهای پرکاربرد در زمینه هیدرولوژی است که فرآیندهای بارش- 
رواناب حوضه‌های آبریز را شبیه‌سازی می‌کند و برای حل گسترده‌ترین مسائل از قبیل هیدرولوژی سیل» تأمين آب 
حوضه رودخانه‌های بزرگ و رواناب حوضه‌های طبیعی يا کوچک شهری طراحی شده است USACE)‏ ۲۰۲۱). 
مدل هیدرولیکی HEC-RAS‏ نيز نسخه ارتقاء‌یافته HEC-2‏ می‌باشد که امکان اجرای > OL‏ ماندگار یک‌بعدی, 


جریان غیرماندگار یک و دو بعدی» محاسبات انتقال رسوب/مرز متحرک و مدل‌سازی دما و کیفیت آب را به كاربر 
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سال یازدهم صحت‌سنجی پهنه و هيدرو GIS‏ سیلاب با استفاده از .... Yo‏ 


می‌دهد (جبلی‌فرد و احمدی, ۱۳۹۵). محققان بسیاری در زمينه استفاده از محصولات مبتنی بر روش‌های دورسنجی 
در مدل‌های هیدرولوژیکی و هیدرولیکی مطالعاتی انجام دادهاند. داسكويتا' و همکاران (۲۰۲۰) در حوضه رودخانه 
ماهانادی " در شمال شرق هند پهنه سیلاب شبیه‌سازی‌شده توسط مدل 1۷11161721 را با پهنه سیلاب استخراج شده از 
تصاویر SAR‏ مقایسه کردند. نتایج» توافق بين پهنه سیلاب مدل‌سازی‌شده و مشاهده‌شده با ضريب 18-0938 را 
نشان دادند. هاریگان " و همکاران (۲۰۲۰) به بررسی داده‌های دبی روزانه بازتحلیل‌شده 010۳۸611۸5 در سطح 
جهان پرداختند. نتایج با مقدار KGE=86%‏ و ميانه ضریب همبستگی پیرسون ۰,۱۱ مناسب ارزیابی شد. پانگ "و 
همکاران (۲۰۲۱) داده‌های بازتحلیل‌شده دبی سیلاب جهانی (GRFR°)‏ را با استفاده از داده‌همای بارش ERAS‏ 
مدل‌های VIC‏ و RAPID‏ در گام‌های زمانی سه‌ساعته و روزانه توسعه دادند. نتايج» عملکرد مناسب مدل‌سازی را در 
همه محدوده‌های = Ob‏ در گام زمانی روزانه و سه‌ساعته در سطح جهان نشان داد. بلابید" و همكاران (۲۰۱۹) در 
منطفه‌ای از اتاوا در کانادا به ارزیابی داده‌های بارش 11۳516۳ و TMPA‏ پرداختند و داده‌های بارش IMERGE‏ را 
به علت دقت SUL‏ به‌عنوان ورودی. به مدل‌های 1۷111 و HEC-RAS‏ وارد کرده و درنهایت پهنه سیلاب حاصل را 
با پهنه سیلاب استخراج‌شده از تصاویر SAR‏ مقایسه نمودند. نتایج» دقت پهنه سیلاب حاصل از داده‌های بارش 
IMERGE‏ را مطلوب نشان دادند. ترینه و ملكنئين" (۲۰۲۱) در حوضه رودخانه ترابانگ" در ویتنام ارزیابی داده‌های 
بارش GSMaP-GNRT6‏ و CMORPH_CRT‏ پرداختند. به علت دقت بالاتر داده‌های بارش CMORPH‏ با ادغام 
مدل‌های MIKE-NAM‏ 211161811 211161821 و MIKE FLOOD‏ پهنه سیلاب آن تهيه شد و بايهنه سیلاب 
استخراج‌شده از تصاویر Sentinel-1A‏ مورد ارزیابی قرار گرفت. نتایج حاکی از دقت مناسب پهنه سیلاب 
مدل‌سازی‌شده نسبت به مشاهداتی بود. در بررسی‌های صورت گرفته در پژوهش‌های پیشین, مطالعات متعددی در 
زمينه اعتبارسنجی داده‌های بارش بازتحلیل‌شده BRAS‏ در سراسر جهان Sl ys)‏ و همکاران ۲۰۱۹ (یوان " و 
همکاران» ۲۰۲۱) و (وروياى'' و همکاران. ۲۰۲۱) و نیز در مناطق مختلف کشور در محدوده‌های استانی 
(عرفانیرحمت‌نیا و همکاران ۱۳۹۷)» حوضه‌ای (گرجی‌زاده و همکاران» ۱۳۹۸) و (عزیزیمبصر و همکاران 
۵) و کل کشور (مانند ایزدی و همکاران ۲۰۲۱) صورت گرفته است؛ اما در زمينه استفاده از این داده‌ها در 
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مدل‌های هیدرولوژیکی به‌ویژه در نرم‌افزار HEC-HMS‏ و بررسی پهنه سیلاب آن در مدل HEC-RAS‏ تحقیقات 
کمی انجام شده اسك بنابراین اهداف پووهش حاضر عبارتند از: ۱) ارزیابی دقت داده‌های بارش یادشده در حوضه 
آبریز کشف‌رود با استفاده از داده‌های بارش ایستگاه‌های باران‌سنجی؛ ۲) بررسی دقت داده‌های بارش BRAS‏ در 
نرم‌افزار HEC-HMS‏ جهت برآورد هیدروگراف جریان در زیرحوضه زشک-برده؛ ۳) بررسی دقت پهنه سیلاب 
برآورد شده در رودخانه زشک با استفاده از دبی پیک هیدروگراف جریان حاصل از داده‌های بارش ERAS‏ در 
راستای اهداف مورد نظر تاکنون مطالعه‌ای در حوضه آبریز کشف‌رود و زیرحوضه زشك-ابرده و رودخانه زشک 
صورت نگرفتهاست. همچنین جهت شبیه‌سازی مشخصات حوضه در نرم‌افزار 110-11015 از ابزارهای موجود در 
منوی GIS‏ استفاده شد كه آوریل ۲۰۲۱ به اين نرم‌افزار افزوده شده است. در واقع اين منو به جای افزونه HEC-‏ 


۹ استفاده می‌شود. 


-Y‏ مواد و روش‌ها 

۱-۲- منطقه مورد مطالعه و داده‌ها 

زیرحوضه زشک- ابرده b)‏ حوضه زشک) در موقعیت جغرافیایی "39 ۵۹74 تا ۱۳ ۵۹۳۱7 طول شرقی و ۱۷ 
۵۲ تا ۱۳ ۲۳ ۳۷۰ عرض شمالی قرار گرفته است (شکل ۱). اين حوضه دارای مساحت حدود ٩۲۲۷,۷‏ 
هکتان محیط ۵۷,4 کیلومتر» ارتفاع متوسط در حدود 774,1 مترء ارتفاع حداقل ۱۶۲۹ و حداکثر ۳۲۲۱ متر و شیب 
متوسط حدود ۳۷,۷ درصد در حوضه آبریز کشف‌رود و از نظر سیاسی در محدوده شهرستان طرقبه- شاندیز واقع 
شده است. طول حوضه ۲۱,۵ کیلومتر بوده و تمامی رواناب اين حوضه از طريق رودخانه زشك-ابرده عبور نموده و 
وارد دشت مشهد می‌شود. بارش متوسط حوضه حدود ۰۶ میلیمتر و بیشترین بارندگی ماهانه به‌ترتیب با ۷۱ و VO‏ 
میلیمتر در ماه‌های فروردین و اردیبهشت اتفاق می‌افتد كه هرکدام حدود ۱۸ و ۱7 درصد از کل بارش را به خود 
اختصاص می‌دهند و ماه‌های خرداد. تیر» مرداد و شهریور نیز خشک‌ترین ماه‌های سال می‌باشند. در پژوهش حاضر 
لایه‌های مورد نياز از اداره كل منابع طبیعی و آبخیزداری استان خراسان رضوی و داده‌های دبی و بارش ایستگاه‌های 
باران‌سنجی حوضه آبریز کشف‌رود از شرکت آب منطقهاى خراسان رضوی دریافت شد که به علت نقص در 
داده‌های ایستگاه‌های حصار-دهبا خرکت. دولت‌آباد خرم‌دره شريف آباد کشف‌رود شغل‌آباد و مزدوران» اين 
ایستگاه‌ها از مراحل انجام پژوهش حذف گردیدند. همچنین داده‌های بازتحلیل‌شده ERAS‏ نيز از پایگاه 
/https://cds.climate.copernicus.eu‏ دریافت گردید. جداول ۱ و ۲ به‌ترتیب حصوصیات فیزیوگرافی حوضه Kij‏ 


و مشخصات ایستگاه‌های باران‌سنجی حوضه آبریز کشف‌رود را نشان می‌دهند. 


سال یازدهم 


صحت‌سنجی پهنه و هيدر وكراف سيلاب با استفاده از .. 


موقميت حوضه زت 


جدول ۱- خصوصیات فیزیوگرافی حوضه زشک (مأخذ: مهندسین مشاور آبخیز گستر شرق) 


نام حوضه 


حوضه آبخیز زشک- 


أبرده 


(km) Lee | مساحت(*16119)‎ 


ov Ad aY YA 


ارتفاع 
حداقل 


(متر) 


۱۶۹ 


ارتفاع 


ارتفاع 


ارتفاع 
میاه 


(متر) 


YYYA 


نمای 
رتفاعی 


(متر) 


110۰ 


طول 


YW 


جدول ۲- موقعیت جغرافیایی ایستگاه‌های باران‌سنجی حوضه آبریز کشف رود JEL)‏ شر کت آب‌منطقه‌ای 


نام ایستگاه 
آبقد ارداک 
آبقد فریزی 
اداره مشهد 
ارداک بند سار وج 
آق دریند 
آل 
امامزاده رادکان 
امامزاده میامی 
اندرخ 
اولنگ اسدی 


خراسان رضوی) 


طول 
04,27 
OA AT‏ 
04,0V‏ 
ERÎ‏ 
VAO‏ 
0417 
04 
7 
011 


0۹۸ 


عرض 
۳۷۲ 
۳,۹ 
۳۳۱ 
AVY‏ 
۳9۹۹ 
ry)‏ 
۳۸ 
۳,۲۳ 
۳,۸ 


۳,۵ 


Migi 
100 


۱۳۹۰ 


YA 


نام ایستگاه 
بلغور 
بهمن جان We‏ 
پل خاتون 
تلغور 
جاغرق 
چکنه عليا 
چناران 
دهانه اخلمد 
زشک خراسان 
سد طرق 
سد کارده 
فریزی 
قدیرآباد 
گلمکان 
كوش بالا 
مارشک 


هندل آباد 


۲-۲- روش پژوهش 


جغرافيا و مخاطرات محیطی 


Migi 
41 


۷۱ 


دقت اندازه گیری بارش برای مدل‌سازی tly,‏ و پیش‌بینی سیلاب‌های شديد و همجنين ساير مخاطرات و 


مورد انتظان داده‌های بارش ERAS‏ يس از دانلود در گام زمانی روزانه و دقت تفکیک ٠ Vx, V‏ برای سال‌های آبی 


1940-7 تا ۱۳۹۸-۹۹ در فرمت ۷600۳ از سایت پادشده در بخش قبلء با استفاده از داده‌های بارش ایستگاه‌های 


باران‌سنجی و شاحص‌های آماری به شرح زيرء با بهره‌گیری از برنامه‌نویسی در محيط ویژوال بیسیک در اکسل 


00 


(Y) 


Xi-1(0i-0)(Pi-P) 


R= 
[BE s(@:-0)? رز‎ eP 


۲ زر‎ Lies (Oi-P) 
NSE’ =1 Si (0;-0) 


1 Trinh و‎ Molkenthin 
2 Nash-Sutcliffe Efficiency (NSE) 


سال یازدهم صحت‌سنجی پهنه و هيدرو گراف سیلاب با استفاده از .... ۳۹ 


w (Y) 
RMSE'= |= )و‎ — 0)? 


)£( م6 یی زووم 

i=1 Vi 

FAR’ = Z (0) 
H 

POD" = = 7 Gj 

15 = FFI K 


که :0 داده‌های مشاهداتی؛ Pi‏ داده‌های شبیه‌سازی‌شده 0 متوسط داده‌های مشاهداتی» H‏ تعداد رویدادهای 
پیش‌بینی‌شده موفقیت‌آمیز توسط یک محصول بارشی معین» F‏ تعداد موارد غیرواقعی که به اشتباه پیش‌بینی شده‌اند و 
M‏ تعداد رویدادهای واقعی كه از دست رفته‌اند. شاحص FAR‏ نسبت هشدار اشتباه که نسبت تعداد رویدادهای 
پیش‌بینی‌شده نادرست به تعداد کل رویدادهای پیش‌بینی‌شده است؛ POD‏ احتمال تشخیص که نسبت تعداد 
رویدادهای پیش‌بینی‌شده درست به كل بارندگی‌های ثبت‌شده توسط ایستگاه‌ها است. TS‏ امتیاز تهدید که اندازه‌گیری 
یکپارچه‌ای از عملکرد کلی را ارائه می‌دهد. محدوده تغييرات اين سه شاخص (۰-۱) است؛ مقدار ايده آل برای 
S25 FAR‏ به صفر و برای TS‏ و POD‏ نزدیک به یک می‌باشد. يس از ارزیابی داده‌های بارش ERAS‏ در مرحله 
بعد لایه مدل رقوم ارتفاعی (DEM)‏ با توان SSH‏ مکانی ۲,۵ امتر ماهواره ALOS-PALSAR‏ برای حوضه زشک 
از وبسایت vertex.daac.asf.alaska.edu‏ دریافت گردید و سپس در نرم‌افزار HEC-HMS‏ فراخوانی شد و با استفاده 
از فرامین موجود در منوی GIS‏ پارامترهای حوضه محاسبه شدند. در ادامه جهت محاسبه میزان تلفات نفوذ از روش 
شماره منحنى SCS‏ استفاده شد. در راستای استفاده ازاین‌روش؛ محاسبات مربوط به درصد نفوذناپذیری و متوسط 
شماره منحنی هر یک از زیرحوضه‌ها به‌ترتیب از طریق clea‏ کاربری اراضی و شماره منحنی به دست آمدند. 
شکل‌های ۲-الف تاج به‌ترتیب لایه‌های گروه هیدرولوژیک خاک و کاربری اراضی و شماره منحنی را نشان 
می‌دهند. به‌منظور تبدیل بارش به رواناب نيز از روش هیدروگراف واحد SCS‏ استفاده شد. يس از معرفی مدل 
حوضه مدل هواشناسی» مشخصه JES‏ و ورود داده‌های بارش و دبى به نرم‌افزار. اقدام به واسنجی مدل شد. بدین 
منظور پس از بررسی وقایع متعدد بارش و رواناب درنهایت دو واقعه مربوط به ٩۱/۱/۱۱‏ و ٩۱/۲/۵‏ جهت واسنجی 


1 Root-Mean-Square Error 
2 False Alarm Ratio 

3 Probability Of Detection 
4 Threat Score 


Yo‏ جغرافیا و مخاطرات محیطی شمارة سوم 


ضريب كارايى نش ساتکلیف (NSE)‏ و درصد 5 استفاده گردیدند. مبنای سنجش شاخص‌های مورداستفاده در 
اين پژوهش, عملكرد شاخصهاى آماری برای مدل‌های حوضه آبریز ارائه شده توسط مورياسى' و همكاران 


N‏ ابه #اربری ا اضر حوصه phd‏ زک ایرد 1 لابه روه رلور یاف حاف حوضه ابر زناف ارده 


راهنما 
گروه هید رولوژیک خاک 
ams‏ 
c a ۱ 1 a pote‏ ل é‏ 4 ۱ 0 
مسکونی i ED Keon eters‏ 


lel, 


agoia شماره‎ 
مقدار‎ 
= a 


ل فا 


۱ 3 a 
مجر نضا‎ ce a ليها‎ Ae `^ 


شکا الف -نقشه گروه هیدرولوژیک خاک حوضه زشک؛ ب- نقشه كاربرى اراضى حوضه زشک؛ ج- نقشه 


شماره منحنی حوضه زشک (دریافتی از اداره کل منابع طبیعی و آبخیزداری استان خراسان رضوی) 


درنهایت DEM a‏ یک متر رودخانه زشک به طول حدود ۳,۳ كيلومتر بالادست ایستگاه زشک از شرکت 
تکاب سبز توس دریافت شد و به همراه دبی پیک سیلاب بدست‌آمده از هیدروگراف سیلاب در مرحله قبلء وارد 
مدل HEC-RAS‏ شده و با استفاده از منوی RAS Mapper‏ خطالقعر رودخانه. سواحل جب و راست» b h>‏ 
مسیر جریان و مقاطع عرضی رسم شدند. در واقع منوی RAS Mapper‏ به جای افزونه HEC-GeoRAS‏ 


1 Moriasi 


سال یازدهم صحت‌سنجی پهنه و هيدرو گراف سیلاب با استفاده از .... ۳۱ 


استفاده می‌شود. روند محاسبه يروفيل سطح آب در حالت جريان ماندگار در مدل هیدرولیکی 3|HEC-RAS‏ یک 
نوع جریان از یکی از دو بازه انتهایی (در جریان فوق بحرانی از مقطع بالادست و در جريان زیربحرانی از مقطع 
avi a22 (A)‏ 


3g و2 + ولا‎ + 29 +h, 


yı +z + 


در این رابطه 1 3 Vo‏ عمق wl‏ در مقطع ابتدایی و انتهایی رودخانه (M)‏ 1 و 2 رقوم کف رودخانه در دو 
مقطع يادشده Vy M)‏ و ج17 سرعت در دو مقطع يادشده m/s)‏ 01 و Ap‏ ضرایب تصحیح انرژی جنبشی در دو 
مقطم يادشده» ع شتاب تقل (۳0/97) و نط افت انرژی (Mn)‏ است. افت انرژی كل نیز از رابطه زیر محاسبه می‌شود: 


av? avi 


2g 2g 


(4) 


hı = رقا‎ +c 


در اين رابطه L‏ فاصله M)‏ و 5۴ شیب اصطکاکی بين دو مقطع و © ضریب افت تنگ‌شدگی و بازشدگی 


] طمن‎ QioptLcn@chtLrop@rob (\ 5 ) 
QlobtQchtQrob 


که در آن gy Liob‏ ترا Len‏ فواصل بين سواحل راست. کانال اصلی و چپ بين هر دو مقطع؛ Qiob‏ و 


Qen « Qrop‏ دبی جریان بين سواحل راست. کانال اصلی و چپ بين هر دو مقطع می‌باشد. درنهایت يس از تعيين 
شرایط مرزی» اجرا و واسنجی مدل» يهنه سيلاب حاصل از دبى ييك بدست‌آمده از هيدر وكراف سیلاب با استفاده از 


داده‌های بارش ERAS‏ (دبی (ERAS‏ سبت به دبى ایستگاه هیدرومتری زشک (دبى مرجع) ارزیابی شد. 
۳- نتایج و بحث 
پژوهش حاضر جهت بررسی پهنه سیلاب حاصل از دبی ERAS‏ نسبت به دبی مشاهداتی انجام شد. بدین 


منظور ابتدا داده‌های بارش 5 در حوضه آبریز کشف‌رود مورد ارزيابى قرار گرفت و پس از واسنجی و 


صحت‌سنجی مدل HEC-HMS‏ با وقايع یادشده در بخش قبلء داده‌های بارش ERAS‏ در موقعیت ایستگاه 


۳۲ 


جغرافیا و مخاطرات محیطی 


زشک مربوط به تاريخ ۹۹/۱/۲۳ وارد HEC-HMS‏ شد و نتایج مربوط به هیدروگراف خروجی أن با نتایج 


مربوط به ایستگاه زشک موردبررسى قرار كرفت. درنهايت يهنه سيلاب حاصل از دبى ERAS‏ نسبت به دبی 


مرجع با استفاده از HEC-RAS‏ و Arc GIS‏ مورد ارزیابی قرار گرفت. در ادامه نتایج مربوط به هر بخش به 


صورت جداگانه آمده است. 


۱-۳- ارزیابی دقت داده‌های باز تحلیل‌شده ERAS‏ در گام زمانی روزانه 


در راستای ارزیابی دقت داده‌های بارش BRAS‏ نسبت به داده‌های بارش زمینی» شاحص‌های آماری مورداستفاده 


در این پژوهش در جدول ۳ آمده است. همانطور که مشاهده می‌شود مقادیر NSE‏ در محدوده مجاز (١-ه)‏ به 


دست آمده است. همچنین RMSE polis‏ و ۸ در محدوده رضایت‌بخشی به دست آمدند. به طور کلی اين محصول 


در تخمین بارش‌های شدید دارای کم‌برآوردی بود که اين نتيجه در مطالعه سینق ' و همکاران (۲۰۲۱) نيز مشاهده 


شده‌بود. شکل ۳ نيز نمودارهای حعبه‌ای شاخص sla‏ محاسبه‌شده در این پژوهش را OLS‏ می‌دهند. 


حدول ۳- جدول polis‏ شاخص‌های آماری جهت ارزیاپی داده‌های بارش ERAS‏ در گام زمانی روزانه 


کت 


TS 


Bias 
RMSE 
POD 


حداقل 


حداکثر 
۷۱ - 
۸44 
I‏ 
۸۲۱ 
YoY‏ 
oV‏ 


۳۲ 


شاخص‌های اماری NSE‏ 00۰۴۸۰6 و Bias‏ 


+ 


FAR 


POD 


NSE 


m a 7 


ميانكين 
۷ - 
FV‏ 


TN 


ميائه 


RMSE مارى‎ 


1 Singh 
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شکل ۳- نمودار جعبه‌ای شاخص‌های TS FAR POD NSE CC‏ و Bias‏ (سمت راست) ونمودار جعبه‌ای 
Co) RMSE‏ جپ) 
۲-۳- ارزیابی دقت داده‌های باز تحلیل‌شده ERAS‏ در گام زمانی ماهانه 
در ارزیابی دقت داده‌های بارش BRAS‏ نسبت به داده‌های بارش زمینی درگام زمانی ماهانه, شاحص‌های آماری 
مورداستفاده در اين پژوهش در جدول ۶ آمده است. شکل ۶ نيز نمودارهای جعبه‌ای شاخص‌های محاسبه‌شده در 
اين پژوهش را OLE‏ می‌دهد. نتایج بدست‌آمده نشان‌دهنده عملکرد نامناسب داده‌های ERAS‏ در ماه‌های گرم و 
خشک (فصل تابستان) بود و در بقيه سال عملکرد مناسب‌تری داشت. به طور کلی polis‏ شاخص‌های ارزیابی به 


مقدار ایده‌آل خود نسبت به گام زمانی روزانه. نزدیک‌تر شده‌اند. 


حدول -٤‏ جدول polis‏ شاخص‌های آماری جهت ارزیابی داده‌های بارش ERAS‏ در گام زمانی ماهانه 


حداقل حداكثر ميانكين میانه 
Mb — y4 Bias‏ = ۳ 10 — 
YW YAYV YAY ۱۷۸۵ RMSE‏ 
“AY Ay AV NÉ POD‏ 
wet wet ۱ 5 FAR‏ 
"M M TS‏ ۷۹ ۷۸ 
AY AY -AY ۸ R‏ 
NSE‏ ۲ وت ۰,00 ۳ ۰,۲۰ 
شاخص‌های آماری «NSE‏ ع۴00۰۴۸8۰6و Bias‏ شاخص آماری RMSE‏ 
یویر ۰ roman‏ تس ۱ ۱ 
0.5 تسس 30 
۰ 10 
کے 
0 
Bias R NSE POD FAR‏ 


شكل -i‏ نمودار جعبه‌ای شاخص‌های TS FAR POD NSE CC‏ و Bias‏ (سمت راست) ونمودار جعبه‌ای 


RMSE‏ (سمت جپ) 


۳-۳- ارزیابی هیدروگراف سیلاب 


Yt‏ جغرافیا و مخاطرات محیطی شمارة سوم 


يس از اعتبارسنجی داده‌های بارش BRAS‏ نرم‌افزار HEC-HMS‏ با استفاده از داده‌های بارش ایستگاه زشک در 
وقايع ۹۱/۱/۱۱ و ۹۱/۲/۱ مورد واسنجی و با استفاده از واقعه مربوط به ۹۹/۱/۲۳ مورد اعتبارسنجی قرار گرفت. 
سپس داده‌های بارش BRAS‏ در گام زمانی روزانه در موقعیت ایستگاه زشک مربوط به واقعه ۹۹/۱/۲۳ وارد HEC-‏ 
5 شد و نتایح Ol‏ با استفاده از شاحص‌های آماری NSE‏ و PBias‏ مورد ارزیابی قرار گرفت. در اين مطالعه 
همانطور که قبلاً بیان شد جهت محاسبه پارامترهای حوضه در نرمافزار 1150-۴1۳15 از منوی GIS‏ شد. شکل ۵ 
شمایی از مدل حوضه ساخته‌شده در نرم‌افزار HEC-HMS‏ با استفاده از منوی GIS‏ را OLS‏ می‌دهد. جدول ۳ نيز 


برخى از پارامترهای محاسبه‌شده توسط نرم‌افزار 1350-11115 را نشان می‌دهد. 


شکل ۵- شمایی از مدل حوضه ساخته شده در HEC-HMS‏ از حوضه زشک 


حدول ۳- حدول برخی از پارامترهای محاسبه شده توسط نرم‌افزار HEC-HMS‏ 


طولانی‌ترین مسیر 
تاش a O‏ شيب ضریب تراکم 
نام dant) | P23‏ مسير جريان 0 تا خروجى طولانی ترين 06 کنیدگی ws‏ 
(km)‏ زیرحوضه‌ها مسير جریان 
(km)‏ 
71 ۱1,۵ 7 كانه ۳,۱ ۳۷ ۳ ۰,۹ ۵۸ 
Ve a ۳,۸ ٤ ye 1۵ 1۹1 72‏ 5۷ 
vY v ۹ £,04 ۱/2 ۷ 11۲ ۱9/۳ 71‏ 
oy £0 0 1,01 ttt ۳ AM ۱۳/۲ 73‏ 
272 ف 1۳ ۸ ۱,0۷ ۳,۷ ۵ ۷ 0 
VE v 3 ۳,۳ \,\0 ۳ RÊ ANA Z3‏ 
"oY ty EAM y «۹ 1,0) qY 74‏ ۰16 
VA 1,07 1,0۳ ۰۷ EY 1,4 7۸‏ ۰5۷ ۳22 


سال يازدهم 


۳۵ .... صحت‌سنجی يهنه و هيدر وگراف سيلاب با استفاده از‎ 
At ity al 6,0 ۳,۹ ۸ 0, Yay Z'5 
VY "oY ۰۳6۵ 0,٤ ۱۳۹۹ «۹ WY 1,۳ Z5 


به منظور مقایسه هیدروگراف‌های مشاهداتی و محاسباتی» شکل‌های 1 تا ۸ در ادامه ارائه شده‌اند. همانطور که 
مشاهده می‌شود مدل در پیش‌بینی بازوی پایین‌رونده هیدروگراف دقت کمتری داشته‌است که اين مسئله به علت 
ورود داده‌های بارش در گام زمانی روزانه (به علت عدم وجود بارانسنج ثبات در حوضه زشک) بود که به تبع آن» 
دقت هیدروگراف محاسباتی نسبت به هیدروگراف مشاهداتی که در گام زمانی ساعتی وارد شده کمتر بوده‌است. 
همچنین مقادیر دبی اوج و زمان رسیدن به دبی اوج نسبت به هیدروگراف مشاهداتی با اختلاف برآورد شده است که 
اين مسئله با نتایج شاحص‌های آماری در جدول ۵ مطابقت دارد. 


شکل -I‏ هیدروگراف واقعه ٩۱/۲/۱‏ 
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جغرافیا و مخاطرات محیطی 


شکل ۸- هیدرو گراف واقعه ۹۹/۱/۲۳ 


حدول ۵-مقادیر حجم و دبى اوج هيدرو گراف‌های مشاهداتى و محاسیاتی بهمراه اختلاف و درصد اختلاف بين 


هیدروگراف 
محاسباتی 
هیدروگراف 
مشاهداتی 
اختلاف 


درصد اختلاف 
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زمان رسیدن 
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QUTO: 
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از بين پارامترهای آماری سنجش خطا نيز در اين پژوهش از ضریب کارایی نش ساتکلیف و درصد Bias‏ استفاده 
شد که نتایج آن در جدول ۱ آمده است. طبق مطالعه موریاسی" و همکاران (۲۰۱۵) در مورد رتبه‌بندی عملکرد 


شاخص هاى ارزیابی برای مدل‌های حوضه آبریزن نتایج مشروحه در جدول ۷ نشان‌دهنده عملكرد مناسب مدل 


هستند زيرا شاخص NSE gl‏ در محدوده‌های بسيارخوب و خوب و PBias‏ در محدوده‌های بسيارخوب» خوباو 


قابل قبول قرار eS‏ در همين جدول مقادير منفى درصد 8135 نشان می‌دهند كه مقادير جريان شبیه‌سازی‌شده 


کمتر از مقادیر مشاهداتی بوده است که اين مسئله همراستا با مطالعه حضریان نژاد و همکاران (۱۳۹۱) می‌باشد. 
به‌طورکلی در همه وقایع» عملکرد مدل رضایت‌بخش بوده است. 


حدول 1- مفادیر شاخص‌های آمارى ضريب كارايى نش ساتكليف و درصد 5 برای وقايع ۱ ۲ ٩۹۱۳۸۱۷"‏ 
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4-۳- ارزیابی هیدروگراف سیلاب با استفاده از داده‌های بارش بازتحلیل شده ERAS‏ 
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سال يازدهم 


صحت‌سنجی پهنه و هيدروكراف سيلاب با استفاده از .... 


در راستای ارزیابی دقت هیدروگراف حاصل از داده‌همای بارش بازتحلیل‌شده ERAS‏ شکل ٩‏ هیدروگراف 
Ob =‏ حاصل از داده‌های دبی ایستگاه هیدرومتری (هیدروگراف مشاهداتی» هیدروگراف جریان حاصل از داده‌های 
بارش ایستگاه زشک (داده‌های بارش زمینی) و هيدر وكراف سیلاب حاصل از داده‌های بارش BRAS‏ را نشان 
می‌دهد. بنابراین شکل» هیدروگراف سیلاب حاصل از داده‌های بارش E845‏ دارای کم‌برآوردی می‌باشد كه علت آن 
کم‌برآورد كردن میزان بارش در موقعیت ایستگاه زشک در محدوده روزهای متناظر با واقعه ۹٩/۱/۲۳‏ توسط ERAS‏ 


بود. عدم تطابق بازوی Gul‏ رونده هیدروگراف‌های داده‌های بارش ایستگاه زشک و BRAS‏ به دلیل عدم برابری 


گام‌های زمانی داده‌های دبى و بارش بود که در بخش قبل توضصیح داده‌شد. جدول ۷مقادیر حجم ودبى اوج 


هيدر وكرافهاى مشاهداتی و محاسباتی به همراه اختلاف و درصد اختلاف بين آن‌ها برای هیدروگراف مشاهداتی و 


هيدر وكراف حاصل از داده‌های بارش ERAS‏ را نشان می‌دهد. 


شکل 4- هیدروگراف سیلاب حاصل از داده‌های ee‏ ایستگاه هید رومتری (منحنی سبز). هيدروكراف سيلاب 
حاصل از داده‌های بارش ایستگاه زشک (منحنی زرد)؛ هیدروگراف سیلاب حاصل از داده‌های بارش ERAS‏ 


(منحنی قرمز) 


حدول -V‏ مقادير حجم و دبى اوج هیدروگراف‌های مشاهداتی و محاسباتی به همراه اختلاف و درصد اختلاف 


بين آن‌ها برای هیدرو گراف مشاهداتی و هیدروگراف حاصل از داده‌های بارش ERAS‏ 


دبی پیک 
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اختلاف 0 اب - 1 24( -VA‏ 1 


درصد اختلاف وت ۷ - AV‏ +- 1 = 


با توجه به مطالعه موریاسی" و همکاران (۲۰۱۵) و شاخحص‌های NSE‏ و درصد Bias‏ در جدول A‏ داده‌های 
بارش بازتحلیل‌شده ERAS‏ با دارابودن صريب PBias NSE‏ در محدوده غير قابل قبول NSES0.5)‏ و 


oak بارشى‎ ‘cme حاصل از داده‌های بارش ایستگاه زشك» كارايى اين‎ ok > نسبت به هيدر وكراف‎ (PBias<+25 


می‌باشد. همچنین با توجه به جدول ۸ به دلیل کم‌برآوردی داده بارش توسط اين محصول. در نرم‌اف زار HEC-HMS‏ 
نيز میزان دبی پیک به مقدار کمتری محاسبه شده است. 


حدول ۸- مقادیر شاخص‌های آماری ضریب کارایی نش ساتکلیف و درصد Bias‏ برای هیدرو گراف مشاهداتی و 


هیدروگراف حاصل از داده‌های بارش ERAS‏ 


شاخص‌های آماري ۳ "مد رو GIS‏ سیلاب مشاهداتی ۳ بارش ERAS‏ 
نش ساتكليف £V E‏ +- 
درصد Bias‏ ۲ - - 


۵-۳- پهنه‌بندی سیلاب 
درنهایت پهنه سیلاب با استفاده از دبی پیک به دست آمده در مرحله قبل به همراه ساير اطلاعات ورودی به 
نرم‌افزار HEC-RAS‏ مورد ارزیابی قرار گرفت. بدین منظورء مدل با استفاده از داده‌های دبی-اشل ایستگاه هیدرومتری 


زشک و تغيير مکرر ضریب زبری مورد واسنجی قرار گرفت. شکل ۱۰ تغییرات مربوط به عمق آب» سرعت جریان 


و تراز سطح آب پس از وا سنجی مدل و اجرای Ol‏ با استفاده از دبی مرجع را نشان می‌دهد. 


شكل ۱۰-تغییرات مربوط به عمق KANT‏ سرعت جریان (ب) و تراز سطح آب (ج) در منوی RAS Mapper‏ 
در نرم‌افزار HEC-RAS‏ 
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در LL‏ پهنه سیلاب به دست آمده از دبی BRAS‏ با پهنه سیلاب به دست آمده از دبی مرجع مورد ارزیابی قرار 
گرفت. شکل ١١‏ پهنه سیلاب حاصل از دبی ERAS‏ و دبی مرجع را نشان می‌دهد. همانطور که مشاهده می‌شود 
وسعت پهنه سیلاب حاصل از دبی ERAS‏ به دليل کم‌برآوردی داده‌های بارش و به دنبال آن کم‌برآورد شدن دبی اوج 
هیدروگراف أن در مدل 11۳0-1115 کمتر از وسعت پهنه سیلاب حاصل از دبی مرجع» به دست آمده است. جدول 


حدول -اختلاف و درصد اختلاف مساحت يهنه سيلاب دبى شبيه سازى شده نسبت به يهنه eo‏ مشاهداتی 
منبع بارشی اختلاف مساحت پهنه سیلاب نسبت به پهنه دبى مشاهداتی درصد اختلاف مساحت پهنه سیلاب نسبت به پهنه دبى مشاهداتی 


-AY -YAYA YY ERA5 


شکل ۱۱-پهنه سیلاب به دست آمده از دبی ERAS‏ و دبی مرجع 


غ- جمع‌بندی 

در پژوهش حاضر به بررسی دقت برآورد پهنه سیلاب حاصل از دبی ERAS‏ نسبت به پهنه سیلاب حاصل از دبی 
مرجع پرداخته شد. بدین منظور ابتدا نیاز به اعتبارسنجی داده‌های بارش اين محصول در حوضه آبریز کشف‌رود بود. 
کم‌برآوردی داده‌های بارش ERAS‏ هم‌راستا با مطالعه ایزدی و همکاران (۲۰۲۱) می‌باشد. همچنین دقت این 
محصول در مدل هیدرولوژیکی HEC-HMS‏ نسبت به هیدروگراف حاصل از داده‌های بارش ایستگاه زشک 
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رضایت‌بخش نبود اما در مطالعه بین ' و همکاران (۲۰۲۰) عملکرد اين محصول قابل قبول بود. پهنه سیلاب دبی 
5 نيز نسبت به پهنه سیلاب دبی مرجع 9 به دست آمد که به دليل کم‌برآوردی داده‌های بارش و به دنبال OF‏ 
کم‌برآورد شدن دبی اوج هیدروگراف آن در مدل HEC-HMS‏ می‌باشد. مطالعات آینده می‌توانند در خصوص 
ارزیابی پهنه سيلاب حاصل از داده‌های بارش محصولات مختلف در نرم‌افزار HEC-RAS‏ باشند. همچنین در جهت 


افزايش دقت اين محصولات. با ترکیب آن‌ها با یکدیگر امجد" و همکاران (۲۰۲۰) می‌توان به اين مهم دست یافت. 
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